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fT^ , Abstract 

o 

. We derive the inclusive transverse spectrum of minijets in nuclear collisions at very high energies, 

^ by assuming that the nuclear ^-matrix factorizes as a product elastic ^-matrices of elementary 

' partonic collisions. Interference effects and, in particular, the contribution of loop diagrams are 

. fully taken into account in the derivation of the spectrum, which is shown to coincide with the 

I result already obtained by superposing the elementary interactions incoherently. A quantitative 

^ • analysis confirms that the deformation induced by multiple coUsions is a large effect at RHIC and 

^ , LHC energies, for transverse momenta ~20 GeV. 

> 

H : 1 Introduction 

The growing experimental activity in heavy-ion coUision has stimulated a renewed interest on the pt 
distributions of secondaries in AB colhsions. From the theoretical point of view a special attractiveness 
presents the production of jets, which, in principle, admits a study in the framework of perturbative 
QCD. In fact theoretical predictions for the inclusive jet production in pp collisions at 1800 GeV agree 
quite well with the experimental data for jet transverse energies in the region 50-h 250 GeV/c[Q. In 
the region of lower jet transverse momenta, the agreement is not as good and phenomenogical param- 
eters need to be introduced|^|^. In nuclear collisions the situation becomes obviously much more 
complicated. There are several mechanisms which influence the form of the jet spectrum and make 
it different from the case of pp. In the nucleus the partonic distributions are changed, the produced 
jets may loose some energy while passing through the thick nuclear matter, jets may accumulate more 
transverse momentum by hard rescatterings. All these different phenomena need to be disentangled 
to gain a comprehensive understanding of the production mechanism. 

In the present paper we focus on the effects of rescattering on the transverse spectrum of minijets in 
nucleus-nucleus collisions. The problem was studied in the past assuming that the total inelastic hard 
cross-section can be obtained from the forward scattering amplitude by substituting the partonic 
(elastic) cross-sections for the forward parton-parton scattering amplitudeQ||]Q, similarly to the 
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analogous case, where the inelastic nucleus-nucleus cross section is expressed in terms of the nucleon- 
nucleon total cross-sections 0, |^. As a consequence, in the approach of Refs. the hard 

nucleus-nucleus cross section is expressed as an incoherent superposition of interactions between the 
partons of the nuclear structures. Possible coherence effects are not included and, in particular, 
the well known problem of loops, arising in the Glauber approach to nucleus-nucleus collisions is 
completely disregarded. 

The purpose of the present paper is to take interferences into account and include all loops in the 
evaluation of rescatterings in the inclusive spectrum of minijet production in nucleus- nucleus collisions. 
In our approach each colliding nucleus will be assumed to consist of a variable number of partons, 
distributed both in impact parameter and rapidity, while the nuclear ^-matrix will be assumed to 
factorize as a product of parton-parton S'-matrices. The inclusive cross-section for minijet production 
will be derived without approximations, apart from the usual disentangling of the longitudinal and 
transverse degrees of freedom and, as we will show, the result will turn out to be identical to the 
expression already obtained in |^ under the incoherent no-loop approximation. 

As an application, in the second part of the paper we study the px distributions in central Pb- 
Pb collisions at nucleon-nucleon cm energies of 200 and 6000 GeV, corresponding to the regimes 
at RHIC and LHC. Since our aim is not to fit the existing data but rather to see the influence of 
rescatterings, we compare our results with the naive picture where each projectile parton is allowed 
to interact with a single target parton only, weighted with the partonic distributions of both nuclei, 
which are taken to be identical to those of the proton. Rescattering effects are obviously maximal 
at relatively small transverse momenta, where the partonic cross-section is relatively large. For this 
reason we concentrate our attention on the production of minijets with transverse momenta up to 10 
20 GeV/c. To make our picture compatible with the existing experimental data in pp collisions at 
this momenta, we shall use an efective partonic cross-sections taken from the lowest order pertubative 
QCD and corrected by phenomenological parameters as in 



2 Glauber theory for hard parton scattering 

We start by introducing the expression for the hard scattering amplitude in terms of partonic inter- 
actions by following the philosophy of the Glauber approach to nucleus- nucleus collisions. Hence we 
assume 1) conservation of partonic longitudinal momenta and impact parameters and 2) factorization 
of the S matrix into the product of elementary partonic S matrices. 

The nucleus state \A) can be represented by a superposition of states with a different number 
n of partons, each characterized by its scaling variable x and impact parameter b. For brevity we 
denote z = {x, b}. We also consider only one partonic flavour, although the generalization to different 
flavours / is trivial. With these notations one has 

/n 
[]d3^,vl/^.„(z,)|n,z,) (1) 

■n i=l 

In the following we shall denote the phase space of the n-parton configuration in nucleus A as 

n 

dTA{n) = l[d^z, (2) 
1=1 

For the nucleus A - nucleus B collision the S matrix has the form 

{A'B'\S\AB) = ^ J dTA{n)dTA{n')dTB{l)dTB{l') 

n.l;n' d' 

*A',„' {^[Wb'^v K)*a.„(z,)*s,;K-)("', r, u]\S\n, z,- 1, u,) (3) 
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Here Uj are the partonic variables of the nucleus B, being I the number of partons in a given configu- 
ration. As stated, our first assumption is that the S matrix is diagonal in the basis {n, Zi; l,Uj}: 

n I 

{n' , z'iJ' ,u'j\S\n, Zi;l,Uj) = 6nn'Su'Y]_S^{zi - zl)Y[S^{ui - u'i)Sni{zi, ...Zn\ui, ...Ul) (4) 

1=1 3=1 

Physically this assumption corresponds taking into account binary collisions only, therefore conserving 
the number of partons in the hard process, which may be justified arguing that in a hard interaction 
parton production is damped by an extra power of the strong coupling constant. We also assume 
that neither the longitudinal momentum nor the impact parameter are changed in a hard collision, 

which arc basic assumptions in the Glauber approach. The range in x and p±^ for which the treatment 
can be applied is therefore restricted to the region where longidudinal momenta are much higher as 
compared to the transferred momenta. Under assumption (4) the S matrix simplifies to 

{A'B'\S\AB) = Y,j dTA{n)dTB{mAAzi)?\'^BM)?Sni{zu -ZnH, ..Ml) (5) 

A second assumption intrinsic in the Glauber approach is that of factorization of the S matrix in 
(5) into a product of elementary S matrices for parton-parton scatterings 

n I 

Snl{zi, ...Zn\ui, ...Ul) = ]J ]J 
i=l j=l 

where Sij is the S matrix for the collision of parton i from nucleus A with parton j from nucleus B: 

Sij = 1 + ia{zi, Uj) = \ + iuij (7) 

a being the parton-parton scattering amplitude. Since we assume that no parton production takes 
place, Sij corresponds to purely elastic scatterings (or to binary reactions in the multi-flavour case). 
Unitarity is therefore satisfied in its simplest form 

{Sij)* Sij = 1 (8) 

(in the multi-flavour case this is changed into a matrix condition in the flavour space) 

To complete the picture we have to introduce the normalization condition for the nuclear wave 
function, which reads 

^ / dTA{n)\^A,n{Zi)\^ = 1 (9) 
n 

As one observes the normalization condition does not limit the parton population of a given nuclear 
configuration, represented by a given term in the sum, it only restricts the total probability of the 
different configurations to be unity. The cross-section introduced in [4] - [6] corresponds to choosing 
the distribution in the number of partons to be Poissonian: 

Wn = J dr(n)|*^,„(^Ol' = e"^"^^ (10) 
If one additionally assumes that the the nuclear partonic wave function factorizes as 

n 

\'^AAZi)\^=Cnl[r{Zi) (11) 
i=l 

one immediately gets 

c„ = ^e-<">, {n)= I d'zTiz) (12) 
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With the forward AB scattering amphtude defined by (5)- (7), we can proceed to determine the 
relevant cross-sections. The total cross-section is of course obtained from the imaginary part of the 
forward {AB 1 5* — 1 1 AB) . However we are not interested in the total cross-section but rather in the total 
hard cross-section with large momenta of observed partons. In the diagrammatic language this cross- 
section corresponds to cutting the forward scattering amplitude through the intermediate partonic 
propagators, with the restriction that the cut partons should have a large transverse momentum 
component. 

One can standardly present all contributions to (5) - (7) in terms of diagrams which show partonic 
interactions via the amplitudes ia.ij. External partons correspond to those present in the colliding 
nuclei. They are assumed to have small transverse momenta and to be distributed in x and b according 
to |^'A(B),nP- AH diagrams can be separated into tree diagrams and diagrams with loops. In the tree 
diagrams all transferred momenta are small (" zero" ) , so that they are expressed through the forward 
scattering amplitudes iaij . Loop diagrams involve integrations over the intermediate parton momenta, 
which may be large. So they cannot be expressed via the forward scattering amplitudes iaij but rather 
involve these amplitudes at large transferred momenta. 

If one restricts oneself to the tree diagrams then the only way to find intermediate partons with 
high transverse momenta is to cut the forward scattering amplitudes iaij converting them into the 
total (elastic) partonic cross-section cr*^'. According to the AGK rules, the sum over all possible cuts 
gives the final hard cross-section in the form {AB\S — l\AB) with iaij substituted by — cfj in (7). 
This latter formula was the starting point in the derivations in Q ||] Q . 

However loops give new contributions to the hard cross-sections, since cutting the loop produces 
intermediate partons with large transverse momenta. The result cannot be expressed in any simple way 
via the partonic cross-sections cr'^', since it involves partonic amplitudes at large momentum transfers. 
Unfortunately it is easy to check that, at least perturbatively, the loop contributions are dominant. 
The reason is that the partonic amplitude a is predominatly real (of order as, whereas its imaginary 
part (and cr'^') is of the order and thus supressed by a factor as). As an explicit example, let us 
compare the different contributions shown in Fig. 1. The external lines of the interacting partons 
supply a factor A^^"^ or B^^^ each, due to the partonic densities in the nuclei. By cutting the single 
scattering diagram a one obtains a contribution of order {AB)^/^a^, while its iteration b gives a 
contribution of order {AB)^/^a^. The contributions of the cuts of the Glauber-like diagrams c and 
d are of orders A^^'^B^^'^a'l and A'^^'^B^^'^a'^. Tree diagrams with 3 interactions, as diagram e in the 
figure, give a contribution of order {AB)^/'^a^. The leading contribution of the loop diagram / is, 
on the contrary of order {AB)^/'^a^ . Since all 4 elementary amplitudes can now be taken uncut (and 
real) in that case, this last contribution is by far dominant, as compared to all connected diagrams, 
and comparable to the iterated lowest order diagram. It can be shown that a similar argument holds 
for any number of interacting partons: the dominant contribution is given by iterations of the sum of 
the lowest order diagrams plus loops. All other diagrams are suppressed by a factor a^. 

As a consequence by substituting ia by — cr'^', which is the approximation used in one 
only selects a limited, sub-leading, set diagrams, the tree diagrams. Notice that the situation is quite 
different as compared to the case of the Glauber picture of AB scattering in terms of nucleons. In 
that case the amplitude is in fact mainly imaginary and, in addition, loops are likely to be suppressed 
by finite formation time arguments. The inclusive hard cross-section, obtained 

in in the 

tree diagram approximation, needs therefore to be re-derived keeping the contribution of loops into 
account. 



5 



3 Jet spectrum: general 



To find the spectrum of emitted jets we liave to fix tlie final state of one of the nuclei in a specific 
manner, in which one of the partons appears with the observed wave function The choice of the 

observed parton is irrelevant, since the total wave function of the nucleus is supposed to be symmetric 
in all the partons, so we choose it to be the first one. Note that choosing the observed parton to 
belong to nucleus A implies studying the spectrum in the direction of nucleus A and conversely for 
nucleus B. In the following we study the former case, so we take 



'bA',n{zi,Z2,...Zn) = ^/n1p{Zl)a^ A' ,n-l{Z2, --Zn) 



(13) 



The second factor represents the state of the remaining n—1 partons, together with the observed one. 
The factor ^/n is to adjust the normalizations of the symmetrized wave functions and ^'^-i. Using 
(13) we get the probability to observe a parton with wave function at a given impact parameter (3 
in the nucleus A-nucleus B scattering as 

^ = ^ ^ /■ dr'A{n')dT'B{l')dTA{n)dTBil) 

^ nl,n',l' A',B'-^ 



C(^l)*Xn(^2,...^n)*B',z(Wj)[5ni(2i|«j) - ll*A,n(2i)*i3,i(Wj) 



(14) 



The sum over all intermediate partonic states, with a given number of partons, is done by using 
closure: 

^ *A,„'(4- ■■■Zn' WA,nizi, --Zn) = ^nn' ^(^'l' ••' ^"1^1' •"' (1^) 
A i=l 

where the 5 function is a symmetrized product of S^[z\ — Zi). We obtain 

^ =^ny" dzidz[i;a{z[)llJa{Zl)dTA{n - IJdTBil)"^ A,n* {z'l, Z2, ...Zn)"^ A,n{zi, Z2, ■■■Zn)]'^ B,l{Uj)\'^ 



S*^l{z[,Z2,...Zn\ui,...,Ui) - 1 Snl{zi,Z2,...Z„\ui,...Ui) - 1 

Notice that by summing over all possible observable states one gets 

Ytpa{z'i)i^a'{zi) = S{z[\zi) 



(16) 



which implies 



E^=E" / dTA{n)drBmA,n\'{Zi)\^BA^^j)\' 
^S*i{z-L,Z2, ...Zn\ui, ...,Ul) - 1 Snl{zi, Z2, ...Zn\ui, ...Ul) - 1 
= ^n / dTA{n)dTB{l)\^A,n\^{Zi)\<S>B,l{Uj)\^ 



nl 



2 - S'*;(zi, Z2, ...Znlui, ...,Ul) - Snl{zi,Z2, ...Z„|ui 



■Ul) 



- <n> (tTb 



(17) 



The result is expected: the integrated inclusive cross-section gives the total cross-section multiplied 
by the average multiplicity, which in our case is obviously < n >. Since < n > is supposed to grow 
with A only linearly (or even less rapidly, if the EMC effect is taken into account) the multiplicity 
does not grow faster then . 
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4 Jet spectrum: non-trivial and geometric contributions 

The product of S matrices in (16) generates 4 terms which we rearrange as 

Snli^'l, Z2, --Znlui, Ul)Snl{zi, Z2, ...Zn\ui, ...Ul) - 1 

- S'nl{z[,Z2,...Zn\Ui,...,Ui)-l - Snl{zi, Z2, ■■■Zn\Ul, ...Ui) - I (18) 

It is instructive to study the meaning of the three terms from the point of view of the intermediate 
physical partonic states, corresponding to cutting the whole diagram with partonic interactions in the 
overall unitarity relation. In the second term only conjugate amplitudes [ia^]* enter, which means 
that this contribution corresponds to cut the interaction diagram from the extreme right, namely to 
the right of all partonic interactions. Conversely in the last term only the amplitudes enter, which 
corresponds to cut the amplitude to the left of all partonic intcgractions. In both cases the momenta 
of the intermediate partons are convoluted with their initial distributions in the colliding nuclei, which 
is very narrow for geometrical reasons, so that also the resulting spectrum cannot be much broader. 

To be more explicit, consider the third term in (18). Take the contribution from a single interaction 
of the first parton at impact parameter hi with some parton in B with its impact parameter c. Since we 
are interested in the distributions of partons from nucleus A, as a function of the transverse momentum 
and at a fixed rapidity, we take 

V'a(^i) = e*f^i(5(.ri-x) (19) 

where x and p are the scaling variable and transverse momentum of the observed parton. Then the 
integration over h and c will take the form 



d'hd'ce'P'^^AAh, ■.■)ia{h - c)|*B,n(c, ..)!' (20) 

where we have suppressed all the variables, with the only exception of bi and c. The typical scale of 
the partonic interaction is much smaller than the nuclear scale, so to a good approximation one can 
rewrite (20) as the product 

J (fhe'P^^^AAbi,...)\^B,n{bi,..)\^ I d^ciaic) (21) 

It is evident now that this expression is different from zero only for values of p of the order of the 
inverse of the nuclear radius, unless one assumes that the nuclear distributions itself contains partons 
with high transverse momenta. 

The non-trivial part of the partonic distribution originated by the hard collisions corresponds 
therefore to cutting the interaction diagram in between hard collisions and thus is totally contained in 
the first term in (18). Remarkably this contribution can be greatly simplified (differently with respect 
to the last two terms in (18)). In fact, using the factorization property (6) we can write 

/ n I 

Sni{zi,Z2,...\ui,...ui) = YlsijYlYls^j (22) 

j=l i=2j=l 

Due to the unitarity relation (8) the second product in (22) is cancelled by its conjugate in the product 
S*S in (18), so that we get 

Snl{Zl,Z2, ...Zn\ui, ...,Ul)Snl{zi,Z2, ...Zn\ui, ...Ul) - 1 = [1 + m(2;i, Wj)]* [1 + «a(^i , Wj)] - 1 (23) 



7 



Putting this expression into (16) we note that (23) does not depend on the partonic variables 
Z2,---Zn of nucleus A. Integrating over Z2,---Zn and summing over n one obtains the p matrix of 
nucleus A as a function of the nuclear partonic degrees of freedom: 



) = pa{zi\z[) 



(24) 



This relation in fact defines the nuclear p-matrix. As a consistency check one may notice that, under 

the factorization assumption (11), at z'^ = zi, (24) reduces to T{zi). 

The contribution from the 1st term in (18) may therefore be written as 

^^(1) „ ; 

-^ = J2 J dzici^;V'a(4)C(^iW(^il4)'^^B(OI*B,iK)l'{n[l + *"(^i'"j)r[l + *"(^i'"i)]-l} 

(25) 

5 The hard spectrum 

To obtain a more explicit form of the spectrum we make use of the factorization ansatz (11) for nucleus 
B: ^ 

l*B,;(^.)l' = ^e-Wnr(«0 (26) 

with 

{l) = j cPuTb{u) (27) 

We also assume for the parton distributions the factorized expression 

Tb{u)=Tb{c)Pb{w) (28) 

where w and c are the fractional momentum and transverse parton coordinate respectively. The 
normalizations are 

j (fbTB{c) = 1, y dwPB{w) = (l) (29) 
Expressing the wave function of the observed parton as in Eq.(19), we obtain 

I I 
JJ (fcjdwjTB{cj)PB{wj){ JJ[1 + ia{z[,Uj)]*[l + ia{zi,Uj)] - l} (30) 

with z = (x, &) and z' = {x,b'). Taking the origin in the center of nucleus B, so that the partonic 
distribution in nucleus A becomes shifted by the overall impact parameter /3, due to factorization, the 
integrations over {wj , Cj ) give the Z-th power of 

J = j (fcdwTB{c)PB{w)[l+ia{x,w;b\-c)]*[l+ia{x,w,bi-c)] (31) 

where the dependence of the scattering amplitude on the transverse distance between the interact- 
ing partons is explicit. The different terms in this integral are treated as in the standard Glauber 
derivation. Let us consider the term with the product of two amplitudes. The distance between the 
interacting partons for hard interactions is very small as compared to the nuclear distances. So using 
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r — bi ~ c as a integration variable we find c = hi — r and we can take the nuclear profile function out 
of the integral at c ~ hi. We obtain 

Tsibi) J (frd'wPB{w)a* {x,w;b[ — bi + r)a{x,w,r) = TB{bi)FB{x,b[ — b) (32) 

where 

(2n 



FB{x,b)^ I d^rdwPBiw)a*ix,w;b + r)a{x,w;r) = j dwPB{w) / 77^1{p)e'^^ (33) 



is the Fourier transform of the transverse momentum distribution in the parton-parton interaction, 
averaged over the parton rapidity distribution of nucleus B. Note that at 5 = (33) gives the total 
cross-section for parton-parton collisions (which coincides with the elastic cross-section), averaged over 
the rapidities of the B-partons 



a'£\x)^ai{x) = j dwPB{w) J d'^r\a{x,w;r)\^ (34) 
The other terms in (31) are evaluated in a similar manner, the result is 

J={1)+ TB{bi)FB{x, b[ -b)- TbFb{x, 0) (35) 
Taking the l-th power of (35) and summing over all values of I one gets 

exp {TB{bi)[FB{xX-hi) - Fb{x,Q)]) 

The term in (30) with unity gives unity after summation over I, so that the cross section is expressed 
as 

{2TTf = / dVb'^e'-P^^'-^'^^ pa{x, 6i - I3\x, b\ - (3) 



d'^Pdyd'^p 

yrB{bl)[FB{x,b[-bt)-FB(x.Q)] _ (^3g) 

Introducing the integration variables r — b[ — bi and b =^ (1/2) (fc^ + 6i) and keeping into account that 
r is a small quantity on the nuclear scale, one obtains 



i'^^f^J^ = /"d26d2re'f>^(a;,6-/3){e^-('')[^«(-''-)-^-(--0)l - l} (37) 



(27r)2_J_ = / d2MVe'^"-T^(6- /3)PA(:^)(e^-(^)[^-(^'^)-^-(^-")l - } (38) 



d-^(3dyd'^p 

Here pa{x, b) = TA{b)PA{x) is just the partonic distribution of nucleus A. Since Fb{x, r) goes to zero 
as r — !■ CX3 the inclusive cross-section (37) contains also a soft component of the spectrum, which in 
the expression is represented by the term proportional to 6^(j>). To remove it we subtract from the 
integrand its value at r = oo. Our final expression for the inclusive cross-section is therefore 
da 
d^Pdyd^p 

The only trace in (38) of the nucleus A is in the term Ta{, b — P)Pa{x), which appears as a weight 
factor for the different contributions from the various parts of nucleus B. For central collisions one 
can separate it as a nearly constant factor, so that the form of the spectrum should be practically 
independent of A and, in particular, it coincides with the spectrum in hadron-nucleus B collisions. 

As mentioned in the Introduction, this formula coincides with the one obtained earlier in [|§§ 
in the tree-diagram approximation. Remarkably, although dominant in the diagrammatic expansion 
of the forward amplitude, as a consequence of unitarity and of the cutting rules, there is no effect of 
loops to the inclusive cross section. 

Notice also that, upon integration over all p, one obtains from (38) the total multiplicity of jets 
from nucleus A multiplied by the cross-section at a given impact parameter (the latter quantity being 
very close to unity for a heavy nucleus). The result coincides with the average number of wounded 
partons of the nucleus A |4|||@§[|o|. From (38) we find 

da 



d^f3dy 



J d^bTA{b-(3)PA{x){l~e"^-'^''^^-^^'''^} (39) 
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6 Numerical results 

The distribution in rapidity and transverse momentum, Eq.(38), is obtained form the elementary 
inclusive cross-section I{u,w, ]p) of two partons colliding with fractional momenta u and w. For the 
gluon- gluon scattering one has ||] 

= ^ [ dy'dudwuG{u)wG{w) ( 1 - —] s(u-x- —] 5 ( w ~ x' ~ —] (40) 
dyd^p 2p^ J \ suw J \ sx' J \ sx J 

with 



P 

—=e", X' — — —e 



aV 

/S ' a/S 



the fractional momenta of the observed and recoiling gluons and G{u), G{w) the distributions of the 
initial gluons, which wc assume to have a small transverse momentum component. Eq. (40) exhibits 
a kinematical constraint between the longitudinal and transverse momenta of the final gluons. In 
particular their scaling variables are somewhat lower than those of the initial gluons. Rigorously 
speaking this goes beyond the Glauber approximation, which, as stated above, assumes the conserva- 
tion of longitudinal momenta. At large s however the problem is felt only close to the limiting value 
~ Pmaa; ~ SUW /A. For p much smaller than its limiting value, at ?/' < one can take G{u) out of 
the integral (40), at a value of momentum fraction corresponding to the momentum of the observed 
gluon, u ^ X. The rest of the integral can be transformed to the integration variable w. In the region 
?/' > one can similarly take G{w) out of the integral and transform the integration to the variable 
u. In this way one obtains the final cross-section in a quasi-Glauber form (see Appendix 1. for the 
details). 

l^='4yGix) I '-^wGiw)(l^^)\(y^-y)] (41) 
dyd-'p p^ I Jwmi„ ^ \ sxw J \ /J 

where 

Wmin = ^(1 + e^y) 

We have checked that the approximations made in reducing (40) to (41) at energies above 200 GeV 
change the spectrum by no more than 10% for p < 40 GeV/c. The two terms in (41) have evidently 
the meaning of the contributions from the projectile and target partons. At y = they are naturally 
equal. To take into account the contribution of quarks, following |H||^[|ll|, we introduce an effective 
parton density P{x) and make the substitution 

xG{x) -> xP{x) = xG{x) + ^ [xQ{x) + xQ{x)] (42) 

Finally, to fit the existing experimental data, we multiply the cross-section by a if -factor. With 
K = 3 and the partonic densities taken from fl^ we obtain a rather satisfactory agreement with the 
experimental data on jet production on protons at 200 GeV [^, as illustrated in Fig. 2. 
The actual expression of the elementary inclusive cross section entering in Eq.(33) is 

2 / „2 " 



I{x,w;p)^{2TrfK^(l-^] e( sxw - p''{l + ey)) (43) 
p^ \ sxw / \ I 



Of course the expression holds only for sufficiently high values of the transferred momenta. So in the 
integral (33) we standardly restrict the integration to p > p^ where po is the infrared cutoff parameter. 
We neglect the change in the initial nuclear parton distributions as compared to the protonic ones 
(the EMC effect) and take Pa{x) — AP(x) where P(x) is the effective partonic distribution in the 
proton (42). In our calculations we use the partonic densities of |12l at the leading order. 
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We restrict ourselves to the production of jets at center rapidity y — 0. To simplify the calculations 
we also consider central collisions {(3 = 0) of identical nuclei and take the nuclei as spheres of radius 
Ra with a constant density, so that 



TA{b) = ^^R\-P (44) 



Then in (38) the integration over b can be done explicitly and we get 

da 



where 



^W = eM4-4 + -)-l-4 (46) 
^ Z'^ z'^ z ) z-^ 



and 

/(.) = -K J .Ve--^(P- ^.^) l^^^ ^-^.PM^)^ [l^f^j (48) 

In (48) we have chosen the scales for as and P{w) to be equal to the transverse momentum squared of 
the recoil parton (i.e. to the square of the transverse momentum transferred to the observed jet). The 
factor 2 in (45) takes into account jets produced by projectile and target partons. For the running 
coupling constant we have taken Aqcd = 0.3 GeV/c with four flavours. 

Some care is needed in determining the scale of the initial partonic distribution P{x). It is related 
to the momentum transferred in the first collision. In the case of a single collision the scale is of order 
p^, the transverse momentum squared of the observed jet. In case of multiple collisions the scale may 
however be substantially smaller. In our calculations we separated the single scattering contribution 
from the general expression (45), where we have taken P{x) at the scale p^. In the remaining part, 
corresponding to double and higher scatterings, we have taken P{x) at the scale p^. 

Even with the above simplifications the calculation of the final inclusive distribution is far from 
trivial since it involves direct and inverse Bessel transforms. We were able to obtain more or less 
stable results only for values of p not higher than ^^20 GeV/c, where, fortunately, we expect hard 
rescatterings to be most relevant. 

We have considered central Pb-Pb coUisons at 200 and 6000 GeV, corresponding to RHIC and 
LHC energies. We have taken po — 2 GeV/c in accordance with the analysis of particle production 
data As for the K-factor we have considered values in the range 1-^3. As mentioned, at 200 
GeV the value A' = 3 is favoured by the proton data The analysis in found that while po 
does not practically change with energy, the K-factor clearly diminishes to values close to unity. 

In Figs. 3 and 4 we present our results in the form of ratios of the distribution (45), with all 
rescatterings included, to the contribution of the single scattering only (optical approximation). The 
latter is obviously just the proton-proton distribution multiplied by a nuclear factor 

7 Conclusions 

Given the very large energies in heavy ion collision at the LHC, a common expectation is that global 
features of the typical event will be within reach of a perturbative QCD approach. It should never- 
theless be emphasized the such a possibility is far from trivial, the capability of making statements on 
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properties of an average inelastic event implying qualitative improvements of present understanding 
of strong interaction dynamics. In a typical heavy ion collision event at the LHC several different 
mechanisms, originated by the complexity of the interactions states, may play important roles giving 
rise to a rather structured interaction, which needs to be brought to light to gain the capability of 
making quantitative statements. 

While an exhaustive description of the whole interaction process, even if limited to the hard 
component only, is still out of reach, various features which are likely to play an important role 
are presently under extensive investigation in the literaturejl^ jl^. The topic which we have 

addressed in the present paper is that of multiparton interactions, where several partons in the initial 
state are linked in a hard collision process, which is a rather natural mechanism to consider in a very 
dense interacting system. One of the reasons of interest is that such a process is able to restore the 
local isotropy in transverse space of the black disk limit of the interaction, which is broken by a two 
parton initiated process at the lowest orders in as [|o| . 

Previous attempts to face this issue were based on a purely probabilistic description of the hard 
component of the multiparton interaction, which has a solid support for disconnected multiparton 
processes, while it has not such a strong basis in the case of a connected multiparton collision. Within 
an incoherent approach, interference effects are in fact neglected altogether, which, in particular, 
implies disregarding all off-diagonal contributions, namely ignoring the well known problem of loops 
in nucleus-nucleus collisions, which give, on the contrary, the dominant contribution to the overall 
interaction amplitude, as discussed in the first part of this paper. To approach the problem we 
have studied the simplest case where only elastic partonic interactions are taken into account and the 
overall nucleus-nucleus S'-matrix is factorized as a product of elementary parton-parton S'-matrices. In 
our approach production processes at the partonic level are therefore completely disregarded. Since 
the complexity of the interaction is described by the Glauber prescription of factorization of the 
S'-matrix, any amplitude for a n-parton interaction process is represented by a convolution of two- 
partons interaction amplitudes, so that the parton-nucleus amplitude is expressed by a series of on 
shell rescattering terms. Notice that, as in the case of the canonical Glauber approach to hadron- 
nucleus collisions, this by no means implies a space-time ordering between different interactions. 
It only implies that a connected n-body interaction process is well approximated by a product of 
two-body interactions, which basically means that, in a dispersive representation of the projectile- 
exchanged gluon amplitude, the pole contribution is dominant. As for the quantitative relevance of 
hard rescatterings, our numerical study confirms that, at transverse momenta of ^ 20 GeV, the effects 
are large both at RHIC and LHC energies. 

The feature which, in our opinion, makes this approach interesting, even if some non-secondary 
features, as energy loss, are not included in the picture of the interaction, is that it allows one to 
obtain, without approximations, the inclusive transverse spectrum of minijets, which is therefore an 
exact consequence of the S'-matrix considered. This feature allows one to argue that the result might 
represent a good starting point for the ambitious program of accomplishing an exhaustive description 
of hard spectra in nucleus-nucleus collisions. 
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8 Appendix 1. Derivation of (41) 

Let us split the integration over y' in (40) in two parts, y' < and y' > 0. Consider the part with 
y' < 0. It corresponds to observed gluons coming from the projectile. Integrating over u we get for 
this part 

ii^) -'41 '41 dwuGiu)raGH (l - ^) % L - - ^) (50) 
with u = X +p'^/{sx'). Further integration over x' gives 

where now 

u = x(l + -^^) (52) 
\ sw — p'^ J 

and Wmin is determined from the condition that x' > pjsfs which gives 

2 

Wmin = —;=-\ (53) 

^^/S XS 

Now we make our crucial approximation to neglect the second term in (53) and put u = x. This 
leads directly to the first term in (41). 

The contribution from the region y' > is transformed in a similar manner. Now we integrate 
over IV to get 

I _ _^ 1 = — ^ / — / dwuG{u)wG{w) (1 — \ S I u- X - — 



u-x--^] (54) 



\dyd'^p J p'^ J x' J \ suw J \ sx' 

with w = x' +p'^/{sx). Further integration over x' gives 

Here 

w=^(l + (56) 

sx \ U — X / 

and Umin is determined from the condition x' < pj \fs to be 

Umin = X -\ p (57) 

Again our approximation consists in neglecting the second term in (56). After that we obtain the 
second term in (41). 

9 Appendix 2. Fits to the existing jet data [7] 

The fits have the form 

da A „h ^gg^ 



dydp (p + Po)" ' GeV/c 
Here follows the table of the parameters with the corresponding values of 
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200 GeV 


500 GeV 


900 GeV 


A 






(^ nn-i- n osW-i-in 
^^i.uuin u.yojii/-riu 


3j 




7 1 1 Q+ 033 




Tin 


4 1 4± n 32 


5 83± 22 


8 50± 28 




3.84 


5.514 


16.4 


A 


(0.91± 0.48)E+8 


(0.175± 0.045)E+8 


(0.379± 0.065)E+7 


a 


6.68± 0.23 


5.383± 0.096 


4.402± 0.068 


Po 


2.0 (fixed) 


2.0 (fixed) 


2.0 (fixed) 




4.54 


10.87 


33.6 



10 Appendix 3. Details of the numerical calculation 

Integration over the angles transforms (45) and (48) in 

da 



(P(3dycPp 



^xAP{x) J rdrJo{pr)(^Z{z) - Z{zo)^ 



and 



where according to (40) and (41) at y — x — p/ ^ and 

2p'' 2p'' 



(59) 



/(r) = 9nK [ p' dp' Up' r)'^^ I —wP{w,p'^) (l - ^\ (60) 



(61) 



^min — — J— 

SX p^Js 

One observes that the two integrations in (60) over w y p' do not decouple, which makes the calculation 
rather complicated. 



11 Figure captions 

1. Diagrams illustrating the interactions of two partons from the projectile with two partons from the 
target. 

2. Jet distributions in pp collisions at 200 GeV at center rapidity. The lower curve is a fit to the 
experimental data [7]. The upper curve shows predictions from Eq. (41) with K = 3.0. 

3. Ratios of the total distributions to the single scattering contributions (optical approximation) 
for central Pb-Pb collisions at 200 GeV and center rapidity calculated from(45) with po = 2 GeV/c. 

4. Ratios of the total distributions to the single scattering contributions (optical approximation) 
for central Pb-Pb collisions at 6000 GeV and center rapidity calculated from(45) dX y — with po = 2 
GeV/c. 
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